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Coulomb gauge Yang- Mills theory within the first order formalism is considered with a view 
to deriving the propagator Dyson-Schwinger equations. The first order formalism is studied with 
special emphasis on the BRS invariance and it is found that there exists two forms of invariance - 
invariance under the standard BRS transform and under a second, non-standard transform. The field 
equations of motion and symmetries are derived explicitly and certain exact relations that simplify 
the formalism are presented. It is shown that the Ward-Takahashi identity arising from invariance 
under the non-standard part of the BRS transform is guaranteed by the functional equations of 
motion. The Feynman rules and the general decomposition of the two-point Green's functions 
are derived. The propagator Dyson-Schwinger equations are derived and certain aspects (energy 
independence of ghost Green's functions and the cancellation of energy divergences) are discussed. 
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INTRODUCTION 



Whilst there is little doubt that Quantum Chromodynamics [QCD] is the theory of the strong interaction, despite 
four decades of intense effort the genuine solution of the confinement puzzle and the hadron spectrum remains elusive. 
That is not to say that no progress has been made - our understanding of QCD is being steadily augmented in many 
ways: for example, with lattice Monte-Carlo techniques [1, 2] and effective theories ([3] and references therein). One 
way to understand the problem of confinement and the hadron spectrum from ab initio principles is to study the 
Dyson-Schwinger equations. These equations are the central equations to the Lagrange formulation of a field theory. 
They are in the continuum and embody all symmetries of the system at hand. 

Dyson-Schwinger studies of QCD in Landau gauge have enjoyed a renaissance in the last ten years. A consistent 
picture of how the important degrees of freedom in the infrared (i.e., those responsible for confinement and the hadron 
""^5 ■ spectrum) stem from the ghost sector of the theory has emerged [4-6] and this has led to increasingly sophisticated 
calculations of QCD properties, culminating recently in hadron observables [7] and possible explanations of confine- 
ment (see the recent review [8] for a discussion of this topic). Landau gauge, in addition to having the appealing 
' property of covariance, has a distinct advantage when searching for practical approximation schemes that allow one 
to extract information about the infrared behaviour of QCD Green's functions, namely, that the ghost-gluon vertex 
remains UV finite to all orders in perturbation theory [9]. For this reason, it has been possible to extract unambiguous 
information about the system [6]. 

Despite the calculational advantages enjoyed by Landau gauge, it is perhaps not the best choice of gauge to study 
the infrared physics of QCD. In this respect, Coulomb gauge is perhaps more advantageous. There exists a natural 
picture (though not a proof) of confinement in Coulomb gauge [10] and phenomenological applications guide the way 
to understanding the spectrum of hadrons, for example in [11, 12] (as indeed they have done in Landau gauge [8, 13]). 
This is not to say that one choice of gauge is better than another - it is crucial to our understanding of the problem 
that more than one gauge is considered: firstly because the physical observables are gauge invariant and it is a test of 
our approximations that the results respect this and secondly because whilst confinement is a gauge invariant reality, 
its mechanism may be manifested differently in different gauges such that we will learn far more by studying the 
different gauges. 

Recently, progress has been made in studying Yang-Mills theory in Coulomb gauge within the Hamiltonian approach 
[14-17]. Here, the advantage is that Gaufi' law can be explicitly resolved (such that, in principle, gauge invariance is 
fully accounted for) and this results in an explicit expression for the static potential between color- charges. In [14-17], 
the Yang-Mills Schroedinger equation was solved variationally for the vacuum state using Gaussian type ansatze for 
the wave-functional. Minimizing the energy density results in a coupled set of Dyson-Schwinger equations which 
have been solved analytically in the infrared [18] and numerically in the entire momentum regime. If the geometric 
structure of the space of gauge orbits reflected by the non-trivial Faddeev-Popov determinant is properly included [16], 
one finds an infrared divergent gluon energy and a linear rising static quark potential - both signals of confinement. 
Furthermore, these confinement properties have been shown to be not dependent on the specific ansatz for the vacuum 
wave- functional but result from the geometric structure of the space of gauge orbits [17]. However, in spite of this 
success, one should bear in mind that an ansatz for the wave-functional is always required and so the approach does 
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not a priori provide a systematic expansion or truncation scheme as, for example, the loop expansion scheme used in 
the common Dyson-Schwinger approach. A study of the Dyson-Schwinger equations in Coulomb gauge will hopefully 
shed some light on the problem. 

Given the appealing properties of Coulomb gauge, it is perhaps surprising that no pure Dyson-Schwinger study 
exists in the literature. However, there is a good reason for this: in Coulomb gauge, closed ghost-loops give rise 
to unregulated divergences - the energy divergence problem. It has been found only relatively recently how one 
may circumvent this problem such that a Dyson-Schwinger study may be attempted [10]. The key lies in using the 
first order formalism. There is not yet a complete proof that the local formulation of Coulomb gauge Yang-Mills 
theory within the first order formalism is renormalisable but significant progress has been made [10, 19]. There is one 
undesirable feature to the first order formalism and this is that the number of fields proliferates. As will be seen in 
this paper, this does have serious implications for the Dyson-Schwinger equations. 

The purpose of the present work is to derive the Dyson-Schwinger equations for Coulomb gauge Yang-Mills theory 
within the first order formalism. These equations will form the basis for an extended program studying QCD in 
Coulomb gauge. The paper is organised as follows. We begin in Section 2 by introducing Yang-Mills theory in 
Coulomb gauge and the first order formalism. In particular, we consider the BRS invariance of the system. Having 
introduced the first order formalism, we then motivate the reasons for considering it (the cancellation of the energy 
divergent sector and the reduction to physical degrees of freedom) in Section 3. Section 4 is then concerned with the 
derivation of the equations of motion and the equations that stem from the BRS invariance. There exists certain 
relationships that give rise to exact statements about the Green's functions that enter the system and these are 
detailed in Section 5. The Feynman rules and the general decomposition of the two-point Green's functions are 
derived and discussed in Sections 6 & 7. In Section 8, the Dyson-Schwinger equations are derived in some detail and 
are discussed. Finally, we summarize and give an outlook of future work in Section 9. 



2. FIRST ORDER FORMALISM AND BRS INVARIANCE 



Throughout this work, we work in Minkowski space and with the following conventions. The metric is g^ v = 
diag(l, —1). Greek letters (fj,, v, . . .) denote Lorentz indices, roman subscripts (i, j, . . .) denote spatial indices and 
superscripts (a, b, . . .) denote color indices. We will sometimes also write configuration space coordinates (x, y, . . .) 
as subscripts where no confusion arises. 

The Yang-Mills action is defined as 



'YM 



'4 P 



where the (antisymmetric) field strength tensor F is given in terms of the gauge field A®: 



F"=d^-d v Al+gr bc A b Al. 



(2.1) 



(2.2) 



In the above, the f abc are the structure constants of the SU(N C ) group whose generators obey [T a ,T b ] = if abc T c . 
The Yang-Mills action is invariant under a local SU(N C ) gauge transform characterised by the parameter 9%,'- 



U x = e^{-i6 a x T a }. 

The field strength tensor can be expressed in terms of the chromo-electric and -magnetic fields (a = A°) 
E a = -d°A a - Va° + , 9 / afcc iV, Bf = , 



VjA% -gf abc A b At 



(2.3) 



(2.4) 



such that Sym = J (E 2 — B 2 )/2. The electric and magnetic terms in the action do not mix under the gauge transform 
which for the gauge fields is written 



- 4 = U x A^Ul - % -{d^U x )Ul 
Given an infinitessimal transform U x = 1 — i9^T a the variation of the gauge field is 



(2.5) 



(2.6) 
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where the covariant derivative in the adjoint representation is given by 

bl c = 5 ac d^ + gf abc Al 



Let us consider the functional integral 



V$ ex-p{iSyM} 



(2.7) 



(2. 



where $ denotes the collection of all fields. Since the action is invariant under gauge transformations, Z is divergent by 
virtue of the zero mode. To overcome this problem we use the Faddeev-Popov technique and introduce a gauge-fixing 
term along with an associated ghost term [20]. Using a Lagrange multiplier field to implement the gauge- fixing, in 
Coulomb gauge (V-A = 0) we can then write 



V$exp{iS Y M + iS/p}, Sf p = I d 



-X a V-A a - c a \7-D ab c b 



(2.9) 



The new term in the action is invariant under the standard BRS transform whereby the infinitessimal parameter 9 a 
is factorised into two Grassmann-valucd components 9 a = c a SX where <5A is the infinitessimal variation (not to be 
confused with the colored Lagrange multiplier field A a ). The BRS transform of the new fields reads 



oc 



1 



-X a SX 
9 

8c a = -lf abc c b c c 5X 
8X a = 0. 



(2-10) 



For reasons that will become clear in the next section (expounded in [10]) we convert to the first order (or phase space) 
formalism by splitting the Yang-Mills action into chromo-electric and -magnetic terms and introducing an auxiliary 
field (7?) via the following identity 



cxpji J d 4 x^E a -E a ^ = J £>7fcxp j« J d 



2 



(2-11) 



Classically, the 7?-field would be the momentum conjugate to A. In order to maintain BRS-invariance, we require that 



J d A x [Sir' 1 - (7?° + £ a ) + Tr a -SE a 



0. 



(2.12) 



Given that the variation of E under the infinitessimal gauge transformation is 8E a = f abc E c 9 h 1 then the general 
solution to Eq. (2.12) is 



abenb 



(1 - a)TT c - aE c 



(2.13) 



where a is some non-colored constant, but which in general could be some function of position x. The 7?-field is split 
into transverse and longitudinal components using the identity 



(2.14) 



const = J V<p8 (V-tt + V 2 

= J V{(j>,T}exp^-i J d i xr a (v-n a + V 2 (f> a 

This constant is gauge invariant and this means that the new fields <ft an d r must transform as 

V 



6cj> a 



(-V 2 ) 



■STT a , Sr a = 0. 



(2.15) 



If we make the change of variables 7? — > 7? — V</> then collecting together all the parts of Z that contain 7?, we can write 



Z 7[ = Jv{n,^,T}ex-pS^i J d 4 x -T a V-n a - ^(jr a - V0 a )-(7r a - V0 a ) - (> - V0 a ) -E a j 



(2.16) 
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which is now invariant under 



5E a 

STT a 



jabc j^CQb 

jabcyb _ a j _ _ 



abc 



(-v 2 ) 



(1 - a) (V - V<£ c ) - aE' 



(2.17) 



We notice that the parts of the transform that are proportional to a arc independent of the rest of the BRS transform 
and can thus be regarded as a separate invariancc. In particular, since it is independent of the Faddeev-Popov 
components, then we may regard it quite generally as a local transform parameterised by 8 a . This new invariance 
stems from the arbitrariness in introducing the 7?-field. If we expand the chromo-electric field into its component form 
then in summary we can write our full functional integral as 



Z = J VQcxpiiSB + iSfp + iSv} 



with 



S B = 


j d 4 x 


2 


1 




Sfp = 


j d 4 x 


-A a V.i* a - 


-c a V-D ab c b 


1 




j d 4 x 


-T a V.7? a " 


- i(7? Q - V0 a )-(7r a - Vcj) a ) + (7? a - V0 a )- (d°A a + D ab a 



and which is invariant under two sets of transforms: the BRS 



5A a 


= -D ac c c SX, 
9 




Sa a 


= --D 0ac c c 5\, 
9 




Sc a 


= -X a SX, 
9 




Sc a 


= -^f abc c b c c SX 7 




Sn a 


= f abc c b Sx(n c - 


V0 C ) + VS<f> 


8(f> a 


fabc J ^ 

" 1 \ (-V 2 ) 




6X a 


= o, 




Sr a 


= o, 





and the new transform, which we denote the a-transform, 

5if a = f abc d b (V - V0 C - d°A - D cd a d ^j + VS(f> a 

S(f> a = f abc 



V 

(-v 2 ) 



• (V - V0 C - d"A c - D cd a c 



d\nb 



(2.18) 



(2.19) 



(2.20) 



(2.21) 



(all other fields being unchanged). It useful for later to denote the combination of fields and differential operators 
occuring in Eq. (2.21) as 



X c = 7? c - V0 C - d"A - D cd cr d . 



(2.22) 
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3. FORMAL REDUCTION TO "PHYSICAL" DEGREES OF FREEDOM 



There are two factors that motivate our use of the first order formalism. The first lies in the ability, albeit formally, 
to reduce the functional integral previously considered (and hence the generating functional) to "physical" degrees of 
freedom [10]. These are the transverse gluon and transverse ir fields which in classical terms would be the configuration 
variables and their momentum conjugates. We keep the term "physical" in quotation marks because it is realised that 
in Yang-Mills theory, the true physical objects would be the color singlet glueballs, their observables being the mass 
spectrum and the decay widths. The second factor concerns the well-known energy divergence problem of Coulomb 
gauge QCD [21-23]. In Coulomb gauge, the Faddeev-Popov operator involves only spatial derivatives and the spatial 
components of the gauge fields, but these fields arc themselves dependent on the spacetime position. This leads to 
the ghost propagator and ghost-gluon vertex being independent of the energy whereas loops involving pure ghost 
components are integrated over both 3-momentum and energy which gives an ill-defined integration. In the usual, 
second order, formulation of the theory these energy divergences do in principle cancel order by order in perturbation 
theory (tested up to two- loops [23]) but this cancellation is difficult to isolate. Within the first order formalism the 
cancellation is made manifest such that the problem of ill-defined integrals can be circumvented. 

Given the functional integral, Eq. (2.18), and the action, Eq. (2.19), we rewrite the Lagrange multiplier terms as 5- 
function constraints and the ghost terms as the original Faddeev-Popov determinant. Since the 5- function constraints 
are now exact we can automatically eliminate any V • A and V • 7? terms in the action. This is clearly at the expense 
of a local formulation and the BRS invariance of the theory is no longer manifest. The functional integral is now 



with 



Z = J 2?$Dct -V-DS A {x-y) S (V- A) S (V-tt) exp {iS} 
= J d 4 x -^B a -B a -\^ a -K a + ^<t> a V 2 <j> a + n a -d a A a +a a (y-D ab (f) b + gp a ^j 



(3.1) 



(3.2) 



where we have defined an effective charge p a — f ade A d ■ ir e . The integral over a can also be written as a S- function 
constraint and is the implementation of the chromo-dynamical equivalent of Gaufi' law giving 



= J 2?$Dct -V-D5 4 (x -y) S (V-/) <5 (v-tt) 8 (-V-D ab <f> b - gp a ^j exp {iS} 



with 



S = J d 4 x 

Let us define the inverse Faddeev-Popov operator M: 



-^B a -B a - ^ a -n a + ^4> a \7 2 4> a + n a -d A a 



ah 



M bc = 5 a 



(3.3) 



(3.4) 



(3.5) 



With this definition we can factorise the GauB law (5-function constraint as 

-l 



(-V-D ab cf) b 



9P 



Det 



-V-D5 A {x-y) ' 5(<j) a - M ab gp b ) . (3.6) 

Crucially, the inverse functional determinant cancels the original Faddeev-Popov determinant, leaving us with 

Z = j V<PS (y-Aj S (V-tt) S (<j) a - M ab gp b ) exp {iS}. (3.7) 

We now use the 5-function constraint to eliminate the (/>-ficld. Recognising the Hcrmitian nature of the inverse 
Faddeev-Popov operator M we can reorder the operators in the action to give us 



Z = j V<P5 (V- A) S (v-n) exp {iS} 



with 



d 4 x 



ls a .B a - ^Tr a -TT a - ^gp b M ba (-\7 2 )M ac gp c + n a -d a A a 



(3.8) 



(3.9) 
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The above action is our desired form, with only transverse A and 7? fields present. All other fields, especially those 
responsible for the Faddeev-Popov determinant (i.e., the certainly unphysical ghosts) have been formally eliminated. 
However, the appearance of the functional (^-functions and the inverse Faddeev-Popov operator M have led to a 
non-local formalism. It is not known how to use forms such as the above in calculational schemes. The issue of 
renormalisability is certainly unclear and one does not have a Ward identity in the usual sense. 

The non-local nature of the above result may not lend itself to calculational devices but does serve as a guide to the 
local formulation. In particular, it is evident that decomposition of degrees of freedom, both physical and unphysical, 
inherent to the first order formalism leads more naturally to the cancellation of the unphysical components in the 
description of physical phenomena than perhaps other choices such as Landau gauge. The task ahead is to identify, 
within the local formulation, how these cancellations arise and to ensure that approximation schemes respect such 
cancellations. For example, the cancellation of the Faddeev-Popov determinant and the appearance of the inverse 
Faddeev-Popov operator should lead to the separation of the physical gluon dynamics contained within the ghost sector 
and the unphysical ghosts themselves, i.e., the unphysical ghost loop of the gluon polarisation should be cancelled 
whilst another loop containing only physical information will take its place. Also, it should be evident that the energy 
divergences associated with ghost loops are explicitly cancelled such that ill-defined integrals do not occur. 



4. FIELD EQUATIONS OF MOTION AND CONTINUOUS SYMMETRIES 



The generating functional of the theory is given by our previously considered functional integral in the presence of 
sources. Explicitly, given the action, Eq. (2.19), we have 



with sources defined by 

S s = J d A x 



Z[J] = J T>$exp{iSB + iSfp + iSn + iSs} 



p a o a + J a -A a + c a r] a + rfc a + K a 4> a + K a -ir a + + &T a 



(4.1) 



(4.2) 



It is useful to introduce a compact notation for the sources and fields and we denote a generic field <& a with source 
J a such that the index a stands for all attributes of the field in question (including its type) such that, for instance, 
we could write 



S s = J a $a (4.3) 

where summation over all discrete indices and integration over all continuous arguments is implicitly understood. 

The field equations of motion are derived from the observation that the integral of a total derivative vanishes 
up to boundary terms. The boundary terms vanish but this is not so trivial in the light of the Gribov problem. 
Perturbatively (expanding around the free-field), there are certainly no boundary terms to be considered, however, 
nonperturbatively the presence of so-called Gribov copies does complicate the picture somewhat. 

In [24], Gribov showed that the Faddeev-Popov technique does not uniquely fix the gauge and showed that even 
after gauge- fixing there are (physically equivalent) gauge configurations related by finite gauge transforms still present. 
It was proposed to restrict the space of gauge field configurations (A) to the so-called Gribov region ft defined by 

ft = |vl : V-A = 0;-V-£> > oj . (4.4) 

ft is a region where the Coulomb gauge condition holds and furthermore, the eigenvalues of the Faddeev-Popov 
operator are all positive. It contains the element A = and is bounded in every direction [25]. However, as explained 
in [26] and references therein, the Gribov region is not entirely free of Gribov copies and one should more correctly 
consider the fundamental modular region A which is defined as the region free of Gribov copies. It turns out though 
that the functional integral is dominated by configurations on the common boundary of A and ft so that, in practise, 
restriction to the Gribov region is sufficient. 

Given that non-trivial boundary conditions are being imposed (i.e., restricting to the Gribov region ft), the question 
of the boundary terms in the derivation of the field equations of motion now becomes extremely relevant. However, by 
definition on the boundary of ft, the Faddeev-Popov determinant vanishes such that the boundary terms are identically 
zero. The form of the field equations of motion is therefore equivalent as if we had extended the integration region to 
the full configuration space [26] . 
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6S S 

6<S>a 



exp {tS + iS s } 



Writing S = Sb + Sf P + Sir we have that 

and so, taking advantage of the linearity in the fields of the source term of the action we have 

J a Z=-J P$ j^jexp^S + iSJ. 



(4.5) 



(4.6) 



We use the convention that all Grassmann-valued derivatives are left-derivatives and so in the above there will be 
an additional minus sign on the left-hand side when a refers to derivatives with respect to either the c-field or the 
77-source. The explicit form of the various field equations of motion are given in Appendix A. 

Continuous transforms, under which the action is invariant, can be regarded as changes of variable and providing 
that the Jacobian is trivial, one is left with an equation relating the variations of the source terms in the action. 
We consider the two invariances derived explicitly in the previous section and that the Jacobian factors are trivial is 
shown in Appendix B. In the case of the BRS transform, Eq. (2.20), we have 



0= I V$—exp{iS + iS s +i5S s } sx=a 



v<s> 



_l a D 0ab c b + }_ja,pab c b 
9 9 



1 



+.f 



abc c b (n c 



v<t> c )- 



K a 



(-V 2 ) 



(n a - V-K a ) 



exp {iS + iS s }. 



(4.7) 



Notice that the infinitessimal variation SX that parameterises the BRS transform is a global quantity, leading to the 
overall integral over x. The equation for the a-transform is: 



/ 



= / V<$>f abc X c x 



K a - 



{K a x -V x -K a x ) 



exp {iS + iS s } 



(4.8) 



where X is given by Eq. (2.22). Constraints imposed by the discrete symmetries of time-reversal and parity will be 
discussed later. 

The above equations of motion and symmetries refer to functional derivatives of the full generating functional. In 
practise we are concerned with connected two-point (propagator) and one-particle irreducible n-point (proper) Green's 
functions since these comprise the least sophisticated common building blocks from which all other amplitudes may 
be constructed. The generating functional of connected Green's functions is W[J] where 



Z[J]=e w W. 

We introduce a bracket notation for functional derivatives of W such that 

SW 



<iJi>- 



81 J 1 



The classical field $ Q is defined as 



1 SZ 

Z SlJ a 



(4.9) 
(4.10) 

(4.11) 



(the classical field is distinct from the quantum fields which are functionally integrated over, but for convenience we 
use the same notation). The generating functional of proper Green's functions is the effective action T, which is a 
function of the classical fields and is defined via a Legendre transform of W: 



r[$] = w[J\ -iJ a $> a . 



(4.12) 



We use the same bracket notation to denote derivatives of T with respect to fields - no confusion arises since we never 
mix derivatives with respect to sources and fields. 

Let us now present the equations of motion in terms of proper functions (from which we will derive the Dyson- 
Schwinger equations). Using the equations of motion listed in Appendix A we have the following equations: 
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• cr-based. This is the functional form of Gaufi' law. 

<uj%> - <ir a x > = V 2 x <t>% + 9f abc A\ x -K c ix - gf abc A b x V lx <i> x + gf abc <iJ b ix iK c ix > 



-9f abc J d 4 yS(x - y)V lx <iJ b iy in x > . 



(4.13) 



Note that we have implicitly used the r equation of motion, Eq. (A. 8), in order to eliminate terms involving 
V-7? in favour of the source £ T . 

• A-based. We write this in such a way as to factorise the functional derivatives and the kinematical factors. The 
equation reads: 

<iAl> = V lx \ a x ^tt? + c4°V„C + [SijVl ~ VixV^] A% 

+g.f abc J d 4 yd 4 z5(y - x)S(z - x) [W iz c b z c c y + n b iz a c y - V iz cj> b z a y ] 

+9f abc J d 4 yd 4 zS(y ~ x)S(z - x) [W iz <i%ir, b z > + <iK b z ip c y > -V lz <in b z ip c y >] 

+9.f abc J d 4 yd 4 z5{y - x)5(z - x) {8 jk W iz + 2%V fc , - S lk V 3V } [<i.J b y i.J c kz > +A b y A c kz ] 

~g 2 f fbc f fde S jk Sli [ S cg S eh( S ab S di + gadgbij + § bg gdh ( § ac gie + § ae § ic^ x 

[<iJ? x iJtiJl> +A% <iJliJl> +A h kx <iJ 9 jx iJl> +A\ X <iJ 9 ]X iJl> +A%A h kx A\ x ] . (4.14) 

• ghost-based. The ghost and the antighost equations provide the same information. The two fields are compli- 
mentary and derivatives must come in pairs if the expression is to survive when sources are set to zero. The 
antighost equation is 

<«£>= -V 2 x c a x - gf abc V lx [<irf x iJ c lx > +c b x A« x ] . (4.15) 



7?-based. 



<<>= V tx r a x - tt? + V«C + d° x A- x + V ix < + gf abc [<ip b x iJ c ix > +° b x A\ x \ 



(4.16) 



0-based. We notice that the interaction terms in the equation of motion for the <^>-field, Eq. (A. 5) are, up to a 
derivative, identical to those of the 7?-based equation, Eq. (A. 6). This arises since — V</> is nothing more than the 
longitudinal part of w and means that there is a redundancy in the formalism that can be exploited to simplify 
proceedings. We can write 



from which it follows that 



(V x -i?« - «S) Z = - 1 V^V 2 x r a x exp {iS + iS s }, 

V x -K a x -n a x = -V£<t£ x >, 
<uj>%> -V x -<ti%> = -V x <. 



(4.17) 

(4.18) 
(4.19) 



• A- and r-based. It will be useful to have these equations written in terms of both connected and proper Green's 
functions. 



e Xx = f x -<if>, <i\ a x >=-v x -i a x , 

& x = V x -<iK£>, <iT a x >=-V x -t a x . 
The BRS transform gives rise to the following equation (the Ward-Takahashi identity): 

= J d 4 x^(d° x p a x ) <vrf x > -f acb p a x [<ip c x ir} b x > + <ip c x ><irf x >] - ^(VixJfJ <ir}> 

-F cb J? x [<iJ? x vrj b x > + <iJ c ix ><vrt>] ~ ~V X «L> ~f abc T x [<%ttf> + <ir&><irf x >] 

9 z 



(4.20) 
(4.21) 



+/ 



abc 



K a — — <K a — V ■ K a \ 



<iK- x irf x > - J d 4 yS(x - y)V ix <iK c x irf y > + <iK c ix ><nf x > - <irf x > V« <m> J . (4.22) 
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We consider for now only the form of the equation relating connected Green's functions. As will be seen in the next 
section, it will not be necessary to consider the equation generated by the invariance under the a-transform. 



5. EXACT RELATIONS FOR GREEN'S FUNCTIONS 



Given the set of 'master' field equations of motion and symmetries, it is pertinent to find out if any of the constraints 
can be combined to give unambiguous information about the eventual Green's functions of the theory. We find that 
such simplifications do in fact exist. 

Let us start by discussing the functional equation generated by a-invariance, Eq. (4.8). It is not necessary here to 
consider functional derivatives of either the generating functional of connected Green's functions (W) or the effective 
action (T) since the derivation applies to the functional integrals directly. From Appendix A, the 7?-based field equation 
of motion, Eq. (A. 6) is 



K? X Z[J] = ~Jv* {V fa 7-° - A7J exp {iS + iS s }, 



Using Eq. (5.1) we can rewrite Eq. (4.8) as 



o = r 



be 



I 



(-vi) 



Kl + V x r! 



exp {tS + iS s }. 



(5.1) 



(5.2) 



Since f abc is antisymmetric and noting Eq. (4.17), the above is now an almost trivial identity. We have thus shown that 
the 7T- and 0-based equations of motion, Eq. (5.1) and Eq. (4.17), guarantee that a-invariance is respected. Conversely, 
approximations to the equations of motion will destroy the symmetry. This is a concrete example of a general feature of 
any physical field theory - the full solutions of the field equations of motion (and the subsequent functional derivatives 
which comprise the Dyson-Schwinger equations) contain all the information given by the symmetry considerations. 
In this case, we have the ambiguity associated with introducing the 7r-field and assigning its properties under the 
BRS transform encoded within the invariance under the a-transform and the field equations of motion are 'aware' of 
this. What is unusual about this, however, is that the equivalence of the field equations of motion and the equations 
generated by invariance under a symmetry is invariably impossible to show (except order by order in perturbation 
theory) - full gauge invariance being the archetypal example. 

Let us now continue the discussion by considering those equations of motion which do not contain interaction terms. 
In the absence of interactions, the solutions to these equations can be written down without difficulty. In terms of 
connected Green's function, the only non-zero functional derivative of the A-equation, Eq. (4.20) is 



(5.3) 



the right-hand side vanishing for all other derivatives. Separating the configuration space arguments and setting our 
conventions for the Fourier transform, we have for a general two-point function (connected or proper) which obeys 
translational invariance: 



<iJ a (y)iJf}(x)>=<iJ a (y - x)iJp(0)>= J dkW a p{k)t 



-ik-{y—x) 



(5.4) 



where dk — d 4 k/(2n) 4 and it is implicitly understood that the relevant prescription to avoid integration over poles 
is present such that the analytic continuation to Euclidean space may be performed. We can immediately write down 
the functional derivatives of <iJf x > using Eq. (4.20): 



<iJ b jy iJ? x > = J dkW b A a Ajl {k)t ]t (k)e 



ik-(y—x) 

*~ " AAjl\' v J"jt\'" / w i 

J dkW b A 3l {k)t 3 S)e- lHy - x \ 



<lK jytj tx^ 



<l£ylJi X > = <iPylJ? x > = <lK v lJ? x >= 



b„ ja 



J 1 „ 7 a ■ 



(5.5) 
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Similarly, for the r-equation, Eq. (4.21), we have 

<iJ b y iK« x > = J dkW^ 3l (k)tdk)e- lHy - x \ 

<iK b ly iKt x > = J dkW% 3l (k)tAk)e- iHy - x \ 

= J dk5 ba ^e-* k <y-*\ 
<ip b y iK? x > = <iK b y iK? x >=<ie Xy iK? x >=0. (5.6) 

We see that, as expected, the propagators involving only the vector fields are transverse and the only other contri- 
butions relate to the Lagrange multiplier fields and are purely kinematical in nature. There is one subtlety to the 
above and that is that whilst the equations for <i£\tJ> and <i£ T iK> are exact in the presence of sources, all other 
equations refer implicitly to the case where the sources are set to zero and the discrete parity symmetry has been 
applied (see later for a more complete discussion). 

In the same fashion, let us consider the A-equation, Eq. (4.20), (and similarly the r-equation, Eq. (4.21)) in terms of 
proper Green's functions. This equation does contain the same information as its counterpart for connected Green's 
functions, but clearly has a different character. The only non-zero functional derivative is 



<iA b „i\ 



>= J dkS ba he- tk < y - x) , (5.7) 



all others vanishing, even in the presence of sources. This applies to all proper n-point functions involving the A-field. 
Similarly, we have the only non-vanishing proper function involving the r-field: 

<«4<>= J dk5 ba k t e- tk <y- x \ (5.8) 

That there are no proper n-point functions involving functional derivatives with respect to the Lagrange multiplier 
fields apart from the two special cases above leads to an important facet concerning the Dyson-Schwingcr equations 
- there will be no self-energy terms involving derivatives with respect to the A- or r-fields since they have no proper 
vertices, despite the fact that the propagators associated with these fields may be non-trivial. 

Next, let us turn to the 0-based equation of motion in terms of connected Green's functions, Eq. (4.18). There are 
only two non-vanishing functional derivatives and we can write down the solutions as before 



<iK b y ig x > = J dkS bai -^e-* k <y-*\ 



I 



<iK b v ie x > = / dks ba ^e- ik <y- xS > 



<<C> = <^C>=<C>=«™C>= 0. (5.9) 

Notice that all of the connected Green's functions involving the T-field are now known and are purely kinematical in 
nature. In terms of proper Green's functions, we consider the 0-based equation of motion, Eq. (4.19). Recognising 
that functional derivatives with respect to the A- and r-fields yield no more information, we can omit them from 
the current discussion. The equation tells us that given a proper Green's function involving n, we can immediately 
construct the corresponding functional derivative with respect to <f>. We can thus conclude that as far as the proper 
Green's functions are concerned, derivatives with repect to the (/>-field are redundant. 

Finally, let us consider the equation derived from the BRS transform in terms of connected Green's functions, 
Eq. (4.22). Since the ghost/ ant ighost fields must come in pairs, we may take the functional derivative of this with 
respect to in d and subsequently set the ghost sources to zero whilst considering only the rest. We get: 



«t>= J d 4 x {- g {dlp%) <irf x inj> -t ch p a x [<ip x ivW z > + <ip><i^t 
-(V ix J? x ) <%rfcrfc> -r b J1 x [<iJZ x irf x ir, d z > + <iJ? x Xitir] d z >] 



9 

j jabc 



K a — Yj± ( K a — V . K a \ 



<iK c lx irf x in d z > - j d A y5{x- y)V ix <in x irf y in d z > + <iKf x Xirf x in d z > - <irf x tr] d > V ix <in c 3 



> 



(5.10) 
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For now, the pertinent information from this identity comes from taking the functional derivative with respect to the 
source £\ and setting all sources to zero. The result is 

«L<L>= o. (5.11) 

We notice that all other functional derivatives lead to non-trivial relations involving interaction terms. This includes 
both the <tpt^\> and the <iki^\> connected Green's functions and we conclude that these functions are not merely 
kinematical factors as one might expect from quantities involving Lagrange multiplier fields. We shall return to this 
topic at a later stage. 



6. FEYNMAN (AND OTHER) RULES 

Whilst it is entirely possible to deduce the complete set of Fcynman rules directly from the action, we shall follow a 
slightly less obvious path here. We derive not only the basic Feynman rules but collect all the tree-level (and incidently 
primitively divergent) quantities that will be of interest. This means that in addition to the tree-level propagators 
(i.e., connected two-point Green's functions) and proper vertices (i.e., proper three- and four-point functions) we 
derive 1 also the proper two-point functions. The reason for this is that (as will be discussed in some detail later) the 
connected and proper two-point functions are not related in the usual way as inverses of one another. The tree-level 
quantities of interest can be easily derived from the respective equations of motion. Indeed, recalling the previous 
section, some are already known exactly and we will not need to discuss them further. 

Before beginning, let us highlight a basic feature of the Fourier transform to momentum space. We know the 
commutation or anti-commutation rules for our fields/sources and this will lead to the simplification that we need 
only consider combinations of fields/sources and let the commutation rules take care of the permutations. However, 
momentum assignments must be uniformly applied and this leads to some non-trivial relations. Consider firstly the 
generic proper two-point function <i$> a (x)i$p(y) > where we have $p(y)$ a (x) = ri$ a (x)$p(y) with r/ = ±1. We 
then have 

<i<f> a (x)i<f>p(y)>= r] <i$p(y)i<$> a (x)> (6.1) 

such that in momentum space 

T a p(k)=r,Tp a (-k). (6.2) 

A similar argument applies for connected two-point functions. The situation for proper three-point functions is slightly 
less complicated since all momenta arc defined as incoming. Indeed, we have (the 5-function expressing momentum 
conservation comes about because of translational invariance) 

<i$ a i<&pi%>= J dk a dkp dk 7 (27r) 4 (5(fc Q + kp + fc 7 )r Q/37 (fc Q , kp, fc 7 ) e -^-*„- t Vx3-**v* 7 ( 6 .3) 
such that, for example 

r / 9a 7 (fc /3 , k a , fey) = r] a/} T a/}l (k a , kp, fey) (6.4) 

where r\ a p refers to the sign incurred when swapping a and 0. 

Let us now consider the connected two-point functions. Setting the coupling to zero in the equations of motion 
(listed in Appendix A) that involve interaction terms gives us the following non-trivial relations (the superscript <>' ^ 
denotes the tree-level quantity) 

Px ~ £tx = — ^x <iK x > ^' 

J?x = -Vi* <^L> (0) +d°x <*tf?><°> ~d° x v lx <!<>("' - [SijVl + v fa v jx ] <j;> (0) , 

Vx = OT7>(°), 

Kl = -V« <C> (0) + > (0) -V fa <<><°> -d° x <iJt x > (0) -V fa < V »><°> . (6.5) 

Clearly, the ghost propagator is distinct from the rest, since the ghost field must appear with its antighost counterpart. 
The treel-level ghost propagator is 

Wfj a \k) = -5 ab ^. (6.6) 
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TABLE I: Tree-level propagators (without color factors) in momentum space. Underlined entries denote exact results. 
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TABLE II: Tree-level proper two-point functions (without color factors) in momentum space. Underlined entries denote exact 
results. Bracketed quantities refer to functions that are fully determined by others. 



The remaining tree-level propagators in momentum space (without the common color factor S ab ) are summarised in 
Table I. Those entries that are underlined are the exact relations considered previously. 

We can repeat the analysis for the tree-level proper two-point functions. The relevant equations are: 

<<><°> - <<><°) = v 2 x r x , 

<iA<i x >W = V lx \ a x -d° x Tr? x + d° x V tx r x +[S l3 V 2 x -V lx V 3X }A« x , 
< t cg>< > = -Vlc a x , 

<*C>(°) - V ix r°-< x + V ix ^ + ^ x + V ix a«. (6.7) 



The ghost proper two-point function is 



f} ab (k) = 5 a \k 2 . (6.8) 



The remaining proper two-point functions are summarised in Table II where the reader is reminded that all proper 
functions involving derivatives with respect to the 0-field can be constructed from the corresponding tt derivative. 

Determining the tree-level vertices (three- and four-point proper Green's functions) follows the same pattern as for 
the two-point functions. They follow by isolating the parts of the equations of motion that have explicit factors of 
the coupling g and functionally differentiating. In momentum space (defining all momenta to be incoming), we have 

r (Q)abc _ _ n rabcc 
1 ncrAij — yJ °ij ' 

T 3Aijk(Pa,Pb,Pc) = ~igf abC [SijiPa ~ Pb)k + Sjk{Pb ~ Pc)t + hi{Pc ~ Pa)j] , 

rfiri? = -w 2 {SijSki [r ce t de - r de r he \ + \j^f^ - r ^f^ + 6il§jk ^^^^ j _ 

4°cm C (^Pc,pa) = -ig.f abc P- a . (6.9) 
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We notice that all the tree-level vertices are independent of the energy. In addition, there is a tree-level vertex 
involving (f> that can be constructed from its counterpart involving 7? and that reads: 

T { 2%<P+P„Pa) = ^f a % = -vf^P*. (6.10) 

This vertex has exactly the same form as the ghost-gluon vertex with the incoming ^-momentum playing the same 
role as the incoming c-momentum. It is worth mentioning that the ghost-gluon, three- and four-gluon vertices are 
identical to the Landau gauge forms except that only the spatial components of the vectors are present. 

Let us now discuss the cancellation of the ghost (energy divergent) sector. In any Feynman diagram containing a 
closed ghost loop, there will be an associated energy divergence. It is a general result that associated with any closed 
loop involving Grassmann- valued fields (ghosts or fcrmions) there will be a factor of (—1). However, Green's functions 
are given by the sum of all possible contributing Feynman diagrams. Since the Feynman rules for W a< p and Y^ aA are 
identical to Wzc and T^cA we will have, for each closed ghost loop, another loop involving scalar fields without the 
factor (—1). Even before performing the loop integration (and regularisation) the integrands of the two diagrams will 
cancel exactly. In this way we see that the energy divergences coming from the ghost sector will be eliminated, as 
expected given that the Faddeev-Popov determinant can formally be cancelled. There is one caveat to this. Whilst 
we have shown that the energy divergences coming from the ghost sector have been eliminated, we have not shown 
that the remaining loops involving scalar fields are free of energy divergences (although a quick glance at the form 
of the Dyson-Schwinger equations later will suffice to see that this is the case at leading order). We propose to look 
further into this in a future publication. 



7. DECOMPOSITION OF TWO-POINT FUNCTIONS 



In order to constrain the possible form of the two-point functions under investigation we can utilise information 
about discrete symmetries. We consider time-reversal and parity and we know that Yang-Mills theory respects both. 
Under time-reversal the generic field & a (x a 7 x) is transformed as follows: 

$ a (x°,x) = 7] a ^ a (-x°,x) (7.1) 

where r\ a = ±1. Since the action, Eq. (2.19), is invariant under time-reversal (it is a pure number) then by considering 
each term in turn, we deduce that 

Va = V\ = 1) Vtt = Vt = V<t> = V<r = -1, i7z = Vc = ±1- (7-2) 

The sources have the same transformation properties as the field. These properties allow us to extract information 
about the energy dependence of Green's functions. For instance, we have that 

Tt ij (k°,k) = -rt ij (-k°,k), (7.3) 

from which one can infer that 

Tf nij (k°, k) = S ab k T Avij (k 2 Q ,k) (7.4) 

(the sign convention is chosen to match the perturbative results). Aside from Ta<p which is unambiguously related to 
T Aw , the only other proper two-point function that carries the external factor k is 

T A b ai (k°,k) = 6 ab k°T Aai (k 2 ,k). (7.5) 

Having extracted the explicit factors of k° in the proper two-point funcitons, the (as yet) unknown functions that 
multiply them are functions of k%. Turning to the propagators, we assign the factor — fc° to Wa v , W a \ and W^x- 
The second discrete symmetry of interest is parity whereby 

<5> a {x\x) = Va $ a (x°,-x) (7.6) 

where again, rj a = ±1. Again the action is invariant and we deduce that 

f7A = % = -l, V<r = V4> = V\ = Vr = 1, r^ = r] c = ±\ (7.7) 

with the sources transforming as the fields. This symmetry is rather more obvious than time-reversal. The physical 
sense is that for every vector field (and with an associated spatial index) we have some explicit vector factor (again 
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TABLE III: General form of propagators in momentum space. The global color factor S ab has been extracted. All unknown 
functions D a p are dimensionless, scalar functions of k% and fc 2 . 



r 






(7 




A 


r 




Uj(k)ik 2 r A A 4- ikikjV A A 




%k Act 


{-ik°ki (r A7T +Ta*)} 


k't 











-ncr 







fcj 


a 


^fc kj~F Act 


— F 


lk Tcra 


$ylk V no | 








4> 


{-ik°kj (r An +T A n)} 




{lk Tttct j- 










A 



















T 





— kj 















TABLE IV: General form of the proper two-point functions in momentum space. The global color factor S ab has been extracted. 
All unknown functions D a p are dimensionless, scalar functions of k% and fc 2 . Bracketed quantities refer to functions that are 
fully determined by others. 



with the associated spatial index). Where the vector fields A and 7? occur in the propagators, we use the transversality 
conditions from earlier to see that the vector-scalar propagators must vanish, except those involving the apropriate 
Lagrange multiplier field. 

What the above tells us is how to construct the most general allowed forms of the two-point functions. The dressing 
functions are scalar functions of the positive, scalar arguments k$ and fc 2 . We summarise the results in Tables III and 
IV. The ghost propagator is written (fe) = —S ab iD c /k 2 . We have nine unknown propagator dressing functions. 

Inlcuding the proper ghost two-point function r^(fc) = 5 ab ik 2 T c we see that there are ten proper two-point dressing 
functions. The extra functions comes about because we have used only the propagator form of the identity Eq. (5.10) 
to eleminate W\\. 

Obviously the propagator and proper two-point dressing functions are related via the Legendre transform. Whereas 
in covariant gauges this relationship is merely an inversion, in our case there is considerably more detail. The 
connection between the connected and proper two-point functions stems from the observation that 

OlJa OlJa OlJ a 

(Recall here that there is an implicit summation over all discrete indices and integration over continuous variables 
labelled by 7.) The ghost two-point functions are somewhat special in that once sources are set to zero, only ghost- 
antighost pairs need be considered. The above relation becomes 

J d 4 z <vrf x iri c z ><ic c z ic b v >= 5 ab S(x - y). (7.9) 
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TABLE V: Possible terms for the equations relating propagator and proper two-point functions stemming from the Legendre 
transform. Entries denote the allowed field types 7 in Eq. (7.8). 



Fourier transforming to momentum space and using the decomposition from above, we get that 

D c {k 2 ,k 2 )Y c {k 2 ,k 2 ) = l (7.10) 

showing that the ghost propagator dressing function is simply the inverse of the ghost proper two-point function. 
Turning to the rest, we are faced with a problem akin to matrix inversion in order to see the connection since 
the sum over all the different possible sources/fields labelled by 7 is non-trivial in the above general formula. The 
decompositions of the two-point function do however mitigate the complexity somewhat. We tabulate the possible 
combinations of terms in Table V. 

We start by considering the top left components of Table V involving only A, 7? and the known functions with A 
and r. After decomposition, we have (suppressing the common argument k) 

klD AlT T A7t - k 2 D AA T AA = k 2 -k 2 , 
k()D A7I T A7T — k D-xftY \k = /cq k , 
D AA T AlT - D A7r T K7r = 0, 

D A kT aa — D n7r T A7r = 0. (7-H) 
We can thus express the propagator functions D in terms of the proper two-point functions T and we have 

1 fc 2 — k 2 1 

D AA = 



(fc 2 - fc 2 ) r AA 



( kg — k J Y A k 

D Av = ^ (7.12) 

Clearly, these expressions can be inverted to give the functions T in terms of the functions D. Next, let us consider 
the central components of Table V involving only a and (f>. We get the following equations: 

^aa^ aa > ^aq}^ na 

Daa^Tia D a( fi (j^nn 4~ ^tttt^) 0, 

Du^Tuu + D^Yna = 0. (7-13) 

The propagator functions in terms of the proper two-point functions are then 



D aa = 



D 1 



r 2 - r (r + r V 



D a A = ^ (7.14) 

r 2 - r (r + r ^ 

X 7T(T L aa \± 7T7T ~ 1 7T7T J 
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There are three more equations that are of interest. These are the a — A, <p — A and A — A entries of Table V and 
they read: 

D aa T Aa - D aii> (r A7r + IV) + D aX = 0, (7.15) 
-D^T Aa - D M (T A7T + T A7T ) + = 0, (7.16) 

T AA - % [D aX T Aa + D^ x (T A7I + T A7I )] = 0. (7.17) 

What these equations tell us is that D a \ and D^x are related to T A7I and T AlJ with all other coefficients being 
determined. is then given as a specific combination and is the 'extra' proper two-point function alluded to 

earlier. However, these functions will not be of any real concern since D a \ and D^x do not enter any loop diagrams 
of the Dyson-Schwinger equations. In effect, T^, T AlJ and form a consistency check on the truncation of the 
Dyson-Schwinger equations since we have that 

- k 2 r - -21 

F AA = -=r- -D ar7 T Aa . + 2D a< f,T A(T (r Aw + T An ) + \T A7r + T Av ) . (7-18) 
k 2 1 J 

It is apparent that unlike covariant gauges, the proper two-point function for the gluon is not necessarily transverse. 

In summary, leaving the problem of the vertices aside, in order to solve the two-point Dyson-Schwinger equations 
we need to calculate seven proper two-point functions: 

Tc: -TAA! TAtT) 1^71-) FtTTT? r<T(T, Tfj-Tr, (7.19) 

which will give us the required propagator functions: 

D c , Daa, D AlT , D V7r , D aa , £> CT 0, (7.20) 

The three proper two-point functions T AA , T A7r and T Aty give a consistency check on any truncation scheme but do 
not directly contribute further. 

8. DERIVATION OF THE PROPAGATOR DYSON-SCHWINGER EQUATIONS 

In this section, we present the explicit derivation of the relevant Dyson-Schwinger equations for proper two-point 
functions. 

A. Ghost Equations 

As will be shown in this subsection, the ghost sector of the theory plays a rather special role. We will begin by 
deriving the ghost Dyson-Schwinger equation (this will serve as a template for the derivation of the other Dyson- 
Schwinger equations). With this it is possible to point out two particular features of the ghost sector: that the 
ghost-gluon vertex is UV finite and that the energy (fc° component) argument of any ghost line is irrelevant, i.e., that 
any proper function involving ghost fields is independent of the ghost energy. 

The derivation of the Dyson-Schwinger equation for the ghost proper two-point function begins with Eq. (4.15). 
Taking the functional derivative with respect to ^c^, using the configuration space definition of the tree-level ghost- 
gluon vertex and omitting terms which will eventually vanish when sources are set to zero, we have 

<ic d w ic? x >= i5 da V 2 J{w - x) + J d 4 yd 4 zT^ c (x iy , z)J^ <zrj b y iJ? z > . (8.1) 

Using partial differentiation we see that 

-^Z < %1 ?y lJ tz>= j d 4 v <iJ9 z irj b y i^><ic e v ici> . (8.2) 

In the above we have again used the fact that when sources are set to zero, the only ghost functions that survive are 
those with pairs of ghost- antighost fields. Since the ghost fields anticommute, we get that 



<.ic x ic w >- 



-i5 ad V 2 J(x -w)+ijd i y d*z d\ T^ c (x, y, z) <iJ^rf y iV><^c d w > (8.3) 
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Taking the partial derivative of Eq. (7.9) with respect to iJf z we have (notice that when using partial derivatives here, 
we must include all possible contributions which for clarity are included explicitly here): 

J d A v <iJ^irf y ivi e v ><ic e v ic d w >^ -i J d 4 vd 4 u <irj b y ir] e v > {<iJliJ f JU ><iA f ]u ic e v ic d w > + <iJ* z iKf u ><VK f ju i? v ici>} . 

(8.4) 

Our ghost Dyson-Schwinger equation in configuration space is thus 
<ic a x ic d w >= -i5 ad X7 2 J(x - w) 

+ j d i yd i zd i vd i uY^i C {x,y,z) <trf y iv> {<iJf z iJ^XiA f ju i^ v ic d w > + <iJ? z iKj u ><mf u ic e v ic d w >} . (8.5) 

We Fourier transform this result to get the Dyson-Schwinger equation for the proper two-point ghost function in 
momentum space: 

T a c d (k) = 6 ad ik 2 

-j{-du) r^f (k, -u,, u> fc)w^H {w c A f Aij (k w )r^.( W) -k, k-u>) + w&^k - w )rf* -k, k-u>)}. 

(8.6) 

With the convention that the self-energy term on the right-hand side has an overall minus sign, we identify (— du; ) 
as the loop integration measure in momentum space. 

With any two-point Dyson-Schwinger equation, it is clear that there are two ordcrings for the functional derivatives 
on the left-hand side. In the same way, there are three orderings for three-point functions and so on. This means that 
there are n different equations for the n-point proper Green's functions, although obviously they all have the same 
solution and must be related in some way. It is therefore instructive to consider also the equation generated by the 
reverse ordering to sec if this will have any consequence. In the ghost case, this means repeating the above analysis 
but starting with the second ghost equation of motion, Eq. (A. 4). The corresponding Dyson-Schwinger equation in 
momentum space is 

Tf (k) = 8 ad ik 2 

-J(-dcj) {w{ c Aij (u>)rg Aj (k, -k + w£yMrg&(fc, -k w, cj)}w b c e (k + u)r^f(k + W) -k, 

(8.7) 

This equation is formally equivalent to Eq. (8.6) but we notice that the ordering of the dressed vertices is different. 
(It is useful to check that the two equations are the same by taking both vertices to be bare such that the equivalence 
is manifest.) 

Notice that one of the vertices that form the loop term(s) must be bare. This arises naturally through the derivation 
above and if one considers a perturbative expansion it is crucial to avoid overcounting of graphs. The choice of which 
vertex is bare is arbitrary and related to the fact that there are n ways of writing the equation for an n-point 
function. Given that for any loop term we can extract a single bare vertex, for any three-point function involving a 
ghost-antighost pair we will have a loop term with the following structure (see also Figure 1): 

/ Swr^.( w ,p, k — p,—k — u)W?(k + u)W{ d aij {u)T^{k + u>, -k, -co). (8.8) 

Now, since the only propagators involving A are transverse (the Wax propagator is disallowed since no proper vertex 
function with A-derivative exists) the loop term must vanish as k — > for finite p[27]. Since the loop term vanishes 
under some finite, kinematical configuration, an UV divergence (which is independent of the kinematical configuration) 
cannot occur and we can can say that this vertex is UV finite. It is tempting to think that such an argument applies 
to the two-point ghost equation, however this is false since whilst the loop term vanishes, so does the k 2 factor that 
multiplies the rest of the equation. 

Let us now show that any Green's function involving a ghost-antighost pair is independent of the ghost and antighost 
energies. The proof of this is perturbative in nature. We notice that both the tree-level ghost propagator and the 
ghost-gluon vertex are independent of the energy. This means that in any one-loop diagram which has at least one 
internal ghost propagator (and hence at least two ghost-gluon vertices) the energy scale associated with the ghost 
propagator is absent. Using energy conservation another energy scale can be eliminated and we choose this to be the 
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FIG. 1: A diagrammatical representation of the T- c/3 (k —p, —k,p) proper vertex dressing. Because of the form of the tree-level 
ghost-gluon vertex and the tranversality of the vector propagator, the dressing function vanishes in the limit k — > 0. Filled 
blobs denote dressed propagators and empty circles denote dressed proper vertex functions. Wavy lines denote proper functions, 
springs deonte connected (propagator) functions and dashed lines denote the ghost propagator. 



antighost energy. At two-loops, we now have the situation whereby the dressed internal ghost propagator is again 
independent of the energy and the dressed ghost-gluon vertex only depends on the gluon energy and so the argument 
can be repeated. This can be applied to all orders in the perturbative expansion which completes the proof. We thus 
have that in particular 

D c {klP) = D c (k 2 ) 
^Acci(ki,k 2 ,h) = T A - Ci (ki,k 2 ,h). (8.9) 



B. (j-based Equation 

Given the discussion in the previous section about which proper two-point functions are relevant, there are only 
two proper two-point functions involving derivatives with respect to a to consider — <ioia> and <tcnir>. Since 
the a-based equation of motion, Eq. (4.13), involves two interaction terms whereas the 7r-based equation, Eq. (4.16), 
has only one, we use the derivatives of Eq. (4.16) to derive the Dyson-Schwinger equation for <iatn> (see next 
subsection). We therefore consider the functional derivative of Eq. (4.13) with respect to ta^, after which the sources 
will be set to zero. We have, again identifying the tree-level vertices, 



<%a d w ia a x >-- 



j d*yd*zr^(z,x,y)^j <iJ b jy iKl> - J d*y d^zT^{z,x,y)-^ oJ» . (8 . 10) 



Using partial differentiation, and with compact notation, 

S 



<ij) ' iKf z >= - <iK c iz iJ a xij] ' iJp><i§ ( )i<§> a io d u} > (8-11) 



(similarly for the second term). This gives the Dyson-Schwinger equation in configuration space: 

<%u d w %ui> = /d 4 !/ d 4 z rg( 2 ,x,j)< l 4 ! j a >< l 4 1 j JJ >< ! $ / , 1 * a ^> 

+ JdSd 4 zT ( £ c A f(z,x,y)<i Kz iJ a ><iJ? y iJp><i$pi$ a iai> (8.12) 
Taking the Fourier transform and tidying-up indices, the Dyson-Schwinger equation in momentum space is thus 
Kt(k) = - J (-^^("-M,-^ 

- j(-dco) rjf( w - k, k, - w )^>)r^ CT >, k co, -k)w£{u> k). (8.13) 

A couple of remarks are in order here. Firstly, there is no bare term on the right-hand side because the action, under 
the first order formalism, is linear in a. Secondly, the implicit summation over the terms labelled by a and (3 means 
that in fact there are eight possible loop terms comprising the self-energy. However, only two of these involves a 
primitively divergent vertex. It is an uncomfortable truth that the formal, non-local delta function constraint arising 
from the linearity of the action in a blossoms into a large set of local self-energy integrals. 
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C. 7r-based equations 

Since the 7r-based equation of motion, Eq. (4.16), contains only a single interaction term, we favor it to calculate the 
<iAnr> and <vki<j> proper two-point functions (as well as <miir>). As discussed previously, we could in principle 
calculate these from the A-based and a-based equations as well in order to check the veracity of any truncations used 
and in fact, this connection may serve useful in elucidating constraints on the form of the truncated vertices used. 
Perturbativcly, all equations will provide the same result at any given order. 

Using the same techniques as in the last subsection, we get the following Dyson-Schwinger equations in momentum 
space: 

KUk) = rt-|(-^)r^^ (8.14) 
K d Mk(k) = -rt% fe -|(-<foOr^^ (8.15) 

r£ jfc (fc) = id ad 5 ik - /(- « W )rg(k, mlm T t«ik("i k ^ -W*M<* - k). (8.16) 

Again, notice that the summation over the allowed types of fields indicated by a and (3 leads to multiple possibilities. 



D. A-based equation 

Using the tree-level forms for the vertices and discarding those terms which will eventually vanish when sources are 
set to zero, it is possible to rewrite Eq. (4.14) as 

<iA1 x > = [SyVl V ix V jx ] A] x + J d 4 yd 4 zr^ a (z,y,x) <irf y z V b z > - J ^dSrg^y,*) <*p» 

-J d^yd^zvfX^y^) -j ^z\v^z, y ,x) <^,^L> 

- J d'y d A z d'w ±rM%j?(w, z, y, x) [3iA b jy <iJ c kz iJf w > +i <iJ b y iJ c kz iJ? w >] • (8.17) 

Functionally differentiating this with respect to A and proceeding as before, noting the following for the four-gluon 
connected vertex 

l KlJjylJ^ljf^ = - <lJ k 2 z ljv><lA e mv t$ i/ i$^><lj fl lj~ f ><ljj y lJx><i®\i$~ f i$ 5 ><lJ S tjf w > 

- <lJk z lJ 1 ><lJ b ]y lJi J ><lA e mv l$ u l§ ll ><lJ^lJ\><l§\l<& 1 l<$>8><lJslJt w > 

- <lJ k 2 Z iJ~f><iJjyiJ\><i®\i$- ( i$5><iJsiJ 'n><lA e mv l§ ^l$„XlJ„ljf w > 

+ <iJl z iJ 1 ><iJ b jy iJ\><iA e mv i<5?xi'& 1 i'&8><iJ8iJt w > (8.18) 
gives the gluon Dyson-Schwinger equation which in momentum space reads: 

rXwW = ^ [k% m hk m ] +J(-du J )r<°£> - k, -w, k)wf(u;)Tg Am (co, k - u , k) 

-J (Su )rg°(, - k, -co, k)W? (w)Tf: Am (u;, k co, -k)W%{w k) 
-J(-du; )Tf b ™{co k, -co, k)W$(u;)Tf: Akm (co, k co, -k)w£ kj (u, - k) 
-\ J ( *<* )r^ (w _ k , -co, k)W c A % t (u;)Tf: Alnm (u;, k - co, -k)W f a b Ank {oo k) 
-I|(-tfu,)(-^)rit,^ 

r{tAno pm (k-CO-V,CO,V,-k) 

+\ J (-dco)Tf A Z? k {k,-k,co,-co)W%\ kl (-uo) 
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FIG. 2: A diagrammatical representation of the coupled system of Dyson-Schwinger equations. Filled blobs denote dressed 
propagators and empty circles denote dressed proper vertex functions. Wavy lines denote proper functions, springs denote 
connected (propagator) functions and dashed lines denote the ghost propagator. Labels indicate the various possible propagator 
and vertex combinations that comprise the self-energy terms. 

+\ J (- )(- dv)T<M$(-v, -u>,v + u>-k, k)W^ ln (v)W^(u)T^ h Xnop (v, w, -v-u)x 

W^ pq (v + oj)T^ Aqrm (v +u,k-v-w, -k)Wi d Arj (u; + v-k). (8.19) 

Again, the occurence of the summation over a, . . . , A leads to many different possible loop terms. 
We present the complete set of Dyson-Schwinger equations in Figure 2. 

9. SUMMARY AND OUTLOOK 

In this work, we have derived the Dyson-Schwinger equations for Coulomb gauge Yang-Mills theory within the first 
order formalism. In discussing the first order formalism it was noted that the standard BRS transform is supplemented 
by a second transform which arises from the ambiguity in setting the gauge transform properties of the n and <fi fields. 
The motivation behind the use of the first order formalism is two-fold: the energy divergent ghost sector can be 
formally eliminated and the system can be formally reduced to physical degrees of freedom, formal here meaning that 
the resulting expressions are non-local and not useful for practical studies. The cancellation of the ghost sector is seen 
within the context of the Dyson-Schwinger equations and the Green's functions stemming from the local action. It 
remains to be seen how the physical degrees of freedom emerge. 

Given that the boundary conditions imposed by considering the Gribov problem and that the Jacobians of both 
the standard BRS transform and its supplemental transform within the first order formalism remain trivial, the field 
equations of motion and the Ward-Takahashi identity have been explicitly derived. The supplemental part of the BRS 
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transform has been shown to be equivalent to the equations of motion at the level of the functional integral and as such 
is more or less trivial. Certain exact (i.e., not containing interaction terms) relations for the Green's functions of the 
theory have been discussed and their solutions presented. These relations serve to simplify the framework considerably. 
The propagators pertaining to vector fields are shown to be transverse, the proper functions involving Lagrange 
multiplier fields reduce to kinematical factors or vanish and the proper functions involving functional derivatives with 
respect to the (/>-field can be explicitly derived from those involving the corresponding 7r-field derivatives. 

The full set of Feynman rules for the system has been derived along with the tree-level proper two-point functions 
and the general form of the two-point functions (connected and proper) has been discussed. The relationship between 
the (connected) propagators and the proper two-point functions, stemming from the Legendre transform, has been 
studied. The resulting equations show that within the first order formalism the dressing functions of the two types 
of two-point Green's functions are non-trivially related to each other. In addition, given that there are no vertices 
involving derivatives with respect to the Lagrange multiplier fields, the set of Dyson-Schwinger equations needed to 
study the two-point functions of the theory is reduced. 

The relevant Dyson-Schwinger equations for the system have been derived in some detail. It is shown how the 
number of self-energy terms is considerably amplified by the introduction of the various fields inherent to the first 
order formalism. The Dyson-Schwinger equations arising from the ghost fields are shown to be independent of the 
ghost energy and the vertices involving the ghost fields are UV finite. 

Despite the complexity of dealing with a non-covariant system with many degrees of freedom, the outlook is positive 
and the rich structure of the Dyson-Schwinger equations is not as intimidating as it might initially appear. The non- 
covariance of the setting means that all the dressing functions that must be calculated are generally functions of 
two variables. However, as seen from the Feynman rules, the energy dependence of the theory stems from the tree- 
level propagators alone and not from vertices (a consequence of the fact that the ony explicit time derivative in the 
action occurs within a kinetic term). The time-dependence of the integral kernels will therefore be significantly less 
complicated than perhaps would otherwise occur. Given the experience in Landau gauge adapting the techniques, 
both analytical and numerical, to solve the Dyson-Schwinger equations in Coulomb gauge seems eminently possible 
though certainly challenging. The results of such a study should provide a better understanding of the issues of 
confinement, and with the inclusion of quarks, the hadron spectrum. 
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APPENDIX A: EXPLICIT FORM OF THE EQUATIONS OF MOTION 



For completeness we list the explicit equations of motion for the various fields represented by Eq. (4.6). 




(A.l) 



J? X Z[J] 



- J V<S> { V„A° - gf bac (V lx c b x )c c x - 8° x « x - V„C) - 9f bac K - Vix4) a% 
+ [SijV x — VtxVjs] Aj x + gf [Aj X \7i X Aj X + 2A^ X SI j X A ix — A ix V jxAj x ~j 

_1„2 ffbc ffde Tzab AC Ae Ad , cad ac Ae Ab , cac Ab Ad Ae , jcae \b Ad At 
J J 1° ^-jx^-jx^-ix ' A jx J± jx Ji ix ' ^-jx^-jx^-ix T" J±j x Ji.j x J± i 



■ix] \exp{iS + iS s }, 



(A.2) 




(A.3) 



(A.4) 



(A.5) 



(A.6) 
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e Xx Z[J] = Jv${v x -A a x }exp{iS + iS s }, 
g x Z[J] = Jv${v x -jr a x }exp{iS + iS s }, 
with X defined in Eq. (2.22). 



(A.7) 
(A.8) 



APPENDIX B: JACOBIAN FACTORS 



The BRS and a-transforms, Eq. (2.20) and Eq. (2.21) respectively, can be regarded as changes of variable in the 
functional integral and that the action is invariant will lead to Ward-Takahashi identities. However, as with all changes 
of integration variable, we must consider the relevant Jacobian factors. These turn out to be trivial, although this is 
not immediately obvious. Starting with the BRS transform, we begin with the observation that 5 A 2 = (since 6 A is 
Grassmann-valued) which means that in the Jacobian determinant only the diagonal elements will contribute since 
all off-diagonal elements are O(SX). Besides the trivial unit terms of the form S ab S(x — y) (with an extra Sij factor for 
vector fields) all diagonal terms have the color structure f abc H c , as can be seen in the form of the transform. Given 
that 



Det(6 ab + f abc H c ) = exp {TrLog(6 ab + f abc H c )} 



(B.l) 



only the first term proportional to f abc survives when the logarithm is expanded (H ~ 5 A). We then see that 
Tr/ abc iJ c = which leaves only the unit term of the exponential. The Jacobian for the BRS transform is thus trivial. 

To see that the Jacobian involved for the a-transform is trivial is slightly more involved. First we note that the 
only non-trivial part of the matrix of variations involves only the 7? and <f> fields, all other rows or columns reducing 
to trivial identity contributions. The sub-matrix of variations for the tt and <f> fields can be written 1 + K where 
K ~ f abc 9 x and is independent of the fields. Since the transform is a change of variables, the functional integral is 
independent of 9% and we can write 



= 



SZ_ 



8=0 



J[6] J V<S> exp {iS + iS s [9}} 



8=0 



With this in mind, it suffices to show that 



6J[6] 



in order for the actual form of the Jacobian to be irrelevant (note that J{9 = 0] = 1). We can write 

J[9] = exp{TrLog(l + K[6])} 



and so 



SJ[0] 



Tr 



8=0 



5K{9] W] SK[0] 



J[9 = 0}. 



(B.2) 



(B.3) 



(B.4) 



(B.5) 



8=0 



Since K is linear in 9, only the first term of the expansion is present when 9 = 0. This term has the color structure 
jabc ^ w j 1 j c i 1 vanishes under the trace operation and so, indeed 



6 J [9] 



= 



(B.6) 



and the Jacobian, although itself not unity is trivial and does not further enter the discussion of the functional integral 
under the transform. 
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